AF IT/GGC/EE/ 78-1 


& 


AN  ANALYSIS  OF  OBSERVED-SYSTEM 
SENSITIVITY  TO  PLANT  PARAMETER 
VARIATIONS 


THESIS 


AF IT/GGC/EE/ 78-1 


Dennis  L.  Hamme 
Capt  USAF 


D D C 

2/ap  nn  n 

JUN  19  1978 

rnsl 

- E 


Approved  for  public  release;  distribution  unlimited 


AFIT/GGC/EE/78-1 


AN  ANALYSIS  OF  OBSERVED-SYSTEM  SENSITIVITY 
TO  PLANT  PARAMETER  VARIATIONS 


THESIS 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science 


ACCESSION  lor 

HUS  ttmto  S*ctiin  91 

C-'C  6un  SmIim  fj 

UNJMIfltlNCFO  q 

JuSllfiMIidH 


BY 

CIl'TNIE  U I ION  IVlUMItlTY  COOtS 


ill:!.  kV,IL  .r  „ or  SPECIAL 

/ 1 

• 

Dennis  L. 


by 

Hamme, 


B . S . E . 


Capt  USAF 

Graduate  Guidance  and  Control 
March  1978 


Approved  for  public  release;  distribution  unlimited 


Preface 

This  topic  was  selected  for  study  as  a result  of  my 
fascination  with  observers.  Although  often  overlooked,  or 
replaced  by  a Kalman  filter,  the  observer  is  a relatively 
simple  solution  to  many  sometimes-complex  problems.  In 
this  report,  a simple  example  plant  was  selected  for  study 
because  I felt  it  would  provide  as  much  insight  to  the 
problem  as  would  a more  comprehensive  model.  Also,  many 
practical  control  systems  are,  in  reality,  no  more  complex 
that  the  two-state  model  investigated.  All  of  the  cases 
worked  were  designed  to  reflect  practical  problems  that 
could  possibly  be  encountered  in  a real-life  situation.  I 
sincerely  hope  that  these  results  will  be  a useful  contribu- 
tion to  the  control  world. 

A special  note  of  appreciation  is  in  order  for  four 
outstanding  persons  who  made  this  project  a little  more 
bearable.  My  lovely  wife,  Peggy,  stood  by  me,  offering 
encouragement  when  everything  seemed  to  go  wrong  and 
rejoiced  with  me  when  obstacles  were  overcome.  My  initial 
thesis  advisor.  Major  Richard  M.  Potter,  helped  immeasur- 
ably in  setting  up  a solid  research  effort  and  getting  me 
started  on  the  problem.  Captain  J.  Gary  Reid  took  over  as 
advisor  well  into  the  study,  picking  up  loose  ends  and 
molding  my  efforts  into  the  final  stretch.  His  patient 

understanding  and  guidance  in  generating  the  computer 
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program  used  for  solving  example  cases  were  of  immeasurable 
assistance  in  accomplishing  the  overall  task.  Finally  a 
note  of  appreciation  is  necessary  for  my  typist,  Phyllis 
Reynolds,  for  her  outstanding  work  and  professional  touch 
that  spared  me  the  normal  hassles  involved  in  typing  a 
thesis.  To  these  four  I am  especially  indebted. 

Dennis  L.  Hamme 
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Abstract 

State-variable  feedback  is  a modern  control  theory 
technique  that  is  employed  in  system  design  to  place 
closed-loop  poles  to  achieve  desired  performance  character- 
istics. Two  problems  associated  with  state-variable  feed- 
back are  physical  inaccessibility  to  plant  states  and  out- 
put sensitivity  to  plant  parameter  variations.  In  this 
report  an  observer  is  employed  to  reconstruct  all  plant 
states.  The  plant-observer  system  is  investigated  with 
respect  to  sensitivity  aspects  of  the  following  three 
areas:  state-variable  representation,  extent  of  pole  place- 
ment, and  observer  dynamics  design.  A comparison  between 
physical,  phase,  and  Jordan  canonical  variables  indicated 
physical  and  phase  variable  representations  yield  identical 
sensitivity  functions;  Jordan  canonical  variables  result 
in  greater  sensitivities  than  physical  and  phase  variables. 
Placing  plant  poles  close  together  creates  a high  system 
sensitivity  to  plant  parameter  variations;  separation  of 
poles  reduces  sensitivity.  As  a nominal  requirement  for 
low  system  sensitivity,  observer  poles  should  be  placed  to 
the  left  of  dominant  plant  poles;  further  reductions  in 
sensitivity  are  achieved  by  placing  observer  poles  to  the 
left  of  all  plant  poles. 
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AN  ANALYSIS  OF  OBSERVED-SYSTEM  SENSITIVITY 
TO  PLANT  PARAMETER  VARIATIONS 

I 

I.  Introduction 


Conventional  control  theory  introduces  the  use  of 
compensators  placed  either  in  cascade  with  or  feedback 
around  basic  plant  elements  to  improve  closed-loop  perform- 
ance by  modifying  the  pole-zero  pattern.  D'Azzo  and  Houpis 
demonstrate  that  the  feedback  of  system  states  through 
appropriate  gains  can  rearrange  the  poles  to  achieve  almost 
any  desired  closed-loop  performance  (Ref  1:410-413).  This 
state-variable  feedback  technique  is  accomplished  by  con- 
structing the  closed-loop  transfer  function,  or  control 
ratio,  in  terms  of  feedback  coefficient,  k^,  for  each  i-th 
system  state.  Using  the  performance  specifications  as  a 
guide,  a desired  control  ratio  is  synthesized  to  achieve 
the  proper  response.  The  final  step  is  to  equate  the  two 
control  ratios  to  determine  the  required  k^.  There  are 
two  major  problems  associated  with  state-variable  feedback 
which  must  be  considered:  all  system  states  may  not  be 
physically  available  for  feeding  back  to  the  input  and 
the  resulting  system  performance  may  be  quite  sensitive  to 
parameter  variations  within  the  plant. 

The  problem  of  inaccessible  states  has  been  investi- 
gated thoroughly  and  several  solutions  proposed.  Minor-loop 
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feedback  of  accessible  states  can  be  employed  to  synthe- 
size inaccessible  states  (Ref  2:517-529).  This  approach, 
however,  requires  exact  knowledge  of  the  plant  parameters 
for  faithful  duplication  of  the  states.  A parameter 
variation  or  inaccuracy  in  determining  plant  parameters 
can  render  this  procedure  impractical.  A second  solution 
is  to  design  an  outer- loop  feedback  compensator  through 
block  diagram  manipulation  (Ref  1:441-444).  This  approach 
requires  utilization  of  differentiator  devices,  which  are 
inherently  noisy,  and  thus  are  not  suitable  for  practical 
applications.  A third  approach  is  to  derive  the  state 
vector  of  a plant  through  modelling  techniques,  such  as 
with  a Kalman  filter  or  an  observer. 

Luenberger  developed  the  mathematical  basis  for  a 
linear  system  that,  when  driven  by  the  inputs  and  outputs 
of  a second  system,  could  reconstruct  the  entire  state 
vector  of  the  second  system  (Ref  3:74-80).  Termed  an 
observer,  this  system  actually  constructs  an  invertible, 
linear  transformation  of  the  state  vector  cr  the  observed 
system  and  has  no  ill  effects  upon  the  control  ratio  other 
than  adding  its  own  poles.  These  are  selected  by  the 
designer  as  a prudent  choice  of  observer  dynamics.  To 
better  understand  this,  the  design  procedure  is  presented 
below.  Given  a linear,  continuous,  time-invariant  system 
of  the  form 

*(t)  * Ax  (t)  + Bu  (t)  (1) 


and  £(t)  = Cx(t)  (2) 

where  x(t)  is  an  nxl  state  vector 
u(t)  is  an  roxl  input  vector 
£(t)  is  an  rxl  output  vector 
A is  an  nxn  transition,  or  dynamics,  matrix 
B is  an  nxm  distribution,  or  control,  matrix 
C is  an  rxn  output  matrix 

an  observer  that  will  reconstruct  a linear  transformation 
of  the  state  vector,  z(t)=Tx(t),  is  given  by 


S(t)  * Dz(t)  + Hx  (t)  + Fu(t) 

(3) 

where 

TA-DT=H 

(4) 

and 

F-TB 

(5) 

The  matrix  H is  selected  to  force  the  observer  dynamics 

matrix,  D,  to  have  some  desired  eigenvalues.  Equations 

(4)  and  (5)  are  then  solved  to  complete  Eq  (3).  The 

matrix  T must  be  invertible  (nonsingular)  to  ensure  the 

*•  _ i 

recovery  of  the  reconstructed  state  vector  x(t)=T  £(t)  . 
There  is  also  a stipulation  that  the  initial  conditions  of 
the  system  and  observer  state  vectors  must  be  equal. 
Luenberger  proved  that  if  the  above  conditions  were  satis- 
fied, the  observer  would  exactly  reproduce  the  selected 
transformation  of  the  system  state  vector  (Ref  3:75). 

In  practical  applications  the  system  output  vector 
rather  than  the  state  vector,  drives  the  observer.  Then 


Eq  (3)  becomes 


4(t)  Dz(t)  + V^ft)  + Fu  (t ) (6) 

or,  by  Eq  (2) , 

z(t)  = Dz(t)  + VCx(t)  + Fu  (t ) (7) 

wherein  Eq  (4)  becomes 

TA-DT=VC  (8) 

This  simplifies  the  calculation  of  D since  only  nr  elements 
of  V must  be  chosen  instead  of  n2  elements  of  H. 

An  observer  is  employed  to  provide  the  inaccessible 
states  for  feedback  to  the  input  to  obtain  the  benefits  of 
state-variable  feedback.  In  fact,  the  entire  reconstructed 
state  vector  can  be  utilized  for  feedback  to  minimize 
instrumentation  costs  by  not  having  to  measure  the  access- 
sible  system  states.  As  mentioned  earlier,  only  the 
observer  poles  are  added  to  the  closed-loop  transfer  func- 
tion and  they  may  be  placed  where  desired  by  the  designer; 
the  plant  poles  are  shifted  independently  by  choice  of 
feedback  coefficients. 

One  significant  question  arises  concerning  the  use 
of  observers  to  enhance  state-variable  feedback:  how  is 
the  control  ratio  affected  by  plant  parameter  variations? 

In  general,  given  any  black  box  with  output  y(t),  the  change 
in  output  due  to  a change  in  an  internal  parameter  is 
expressed  by  a first-order  Taylor  series  approximation  as 
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y(t;p0+Ap)  : y(t;pQ)  + 8y-'-apf- 


Ap  (9) 


where  pQ  is  the  nominal  parameter  value  and  Ap  is  the 
variation.  If  the  system  is  linear,  the  steady-state  out- 
put can  be  considered  with  no  loss  of  generality.  Rewriting 
Eq  (9)  as 


yss(t;p0+Ap) 


y (t;p  ) 

■'ss  ^o 


3yss(t;p) 

ap 


(10) 


it  is  seen  that  the  change  in  output  is  approximately  equal 
to  the  "sensitivity"  of  y(t)  to  p scaled  by  the  variation 
in  p.  Since  the  sign  and  magnitude  of  Ap  are  not  known, 
the  area  of  concern  lies  in  the  magnitude  of  the  sensitivity 
function.  For  greater  understanding  of  how  the  system  is 
affected  by  Ap,  a frequency  response  of  the  sensitivity 
function  would  reveal  any  particular  frequencies  of  opera- 
tion for  which  the  system  would  be  especially  susceptible 
to  Ap. 

Consider,  for  example,  the  linear  regulator  in 
Fig.  1.  The  output  of  the  plant  is  fed  back  to  nullify 
disturbance  inputs,  d (t) . If 
the  output  is  particularly 
sensitive  to  a certain  band  of 
frequencies  within  the  opera- 
tional bandwidth,  the  feedback 
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signal  could  itself  introduce  disturbance  inputs  for  some 
parameter  variation  within  the  plant. 

Another  example  is  in  the  control  system  in  Fig.  2. 


Figure  2.  Feedback  Control  System 


Utilizing  state-variable  feedback,  an  input  r(t)  causes  a 
desired  steady-state  output  y(t).  The  feedback  coefficient 
matrix,  K,  is  selected  to  provide  a particular  response  and 
the  observer  constructs  the  required  state  vector.  Under 
a parameter  variation  the  plant  output  may  vary,  causing 
errors  in  the  observer  which,  in  turn,  generates  the  wrong 
feedback  signal  that  creates  more  errors — and  the  cycle 
continues.  Or,  quite  possibly,  the  plant  parameters  are 
not  known  exactly  and  the  observer  design  is  based  upon 
"best  guess"  parameters.  The  reconstructed  state  vector 
will  contain  an  error  such  that  the  feedback  signal  does 
not  provide  the  desired  system  response.  In  both  cases 
the  sensitivity  function  provides  insight  into  the  severity 
of  output  deviations  due  to  parameter  variations,  A fre- 
quency response  of  the  sensitivity  function  can  assist  in 
determining  whether  or  not  the  system  can  perform 
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satisfactorily  within  its  operating  bandwidth. 

The  objective  of  this  thesis  is  to  compare  the 
sensitivity  functions  for  direct  state  feedback  and 
observer-implemented  feedback  systems  to  determine  if  the 
employment  of  observers  is  acceptable  in  light  of  plant 
parameter  variations. 


II.  Development  of  the  Study 


There  are  three  major  areas  to  be  investigated 
for  their  effects  upon  system  sensitivity  to  parameter 
variations:  state-variable  representation,  pole  placement, 
and  observer  design.  These  topics  will  be  studied  both 
individually  and  in  combinations  to  test  for  possible 
interdependence . 

It  has  been  demonstrated  that  the  choice  of  state- 
variable  representation  can  have  a rather  drastic  effect 
upon  eigenvalue  sensitivity.  A small  change  in  a plant 
parameter  can  cause  the  eigenvalues  to  vary  by  many  orders 
of  magnitude  or  very  little  at  all  (Refs  4:85-93?  5:263- 
266) . The  concept  of  eigenvalue  sensitivity  centers 
around  simulation  of  a linear,  time-invariant  system  on  a 
computer.  A high  sensitivity  indicates  an  eigenvalue  of 
the  system  may  vary  greatly  from  its  nominal  value  in  the 
presence  of  elemental  variations  or  inaccuracies  in  model- 
ling. The  resulting  errors  in  system  dynamics  detract 
from  the  benefits  of  modelling.  The  correlation  between 
eigenvalue  sensitivity  and  closed-loop  system  sensitivity 
to  parameter  variations  is  not  well-defined?  therefore, 
a comparison  will  be  made  between  three  common  types  of 
state  variables  to  determine  whether  the  type  of  state- 
variable  representation  affects  system  sensitivity. 

Physical  variables  are  selected  because  they  depict  actual. 
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measurable  quantities  that  can  be  utilized  for  feedback. 
Phase  canonical,  or  simply  phase,  and  Jordan  canonical 
variables  are  purely  mathematical  quantities  that  have 
practical  applications  in  control  system  analysis  and 
design.  Singer  shows  that  these  two  forms  can  have 
radically  different  eigenvalue  sensitivities  (Ref  4:88) 
and,  thus,  are  logical  candidates  for  study. 

The  second  area  of  interest  is  the  pole  placement 
that  results  from  state-variable  feedback.  It  seems 
possible  that  the  distance  poles  are  shifted  and/or  the 
pattern  desired  (e.g.,  double  roots)  may  affect  system 
sensitivity  because  of  the  range  of  possible  feedback 
gains.  An  initial  guess  is  that  larger  feedback  gains 
would  tend  to  decrease  system  sensitivity  due  to  increased 
feedback  effects,  except  that  errors  in  the  reconstructed 
state  vector  would  be  amplified.  The  converse  seems 
logical  for  small  feedback  gains.  This  facet  of  the 
problem  will  be  studied  to  determine  if,  indeed,  there 
exists  any  correlation  between  pole  placement  and  sensi- 
tivity. 

There  have  been  many  articles  published  on  observer 
theory  (see,  for  example  Refs  7,  8,  and  9)  but  none  of  them 
specifically  address  general  guidance  in  observer  pole 
selection.  Luenberger  states  only  that  the  observer 
dynamics  should  have  a quicker  response  than  the  dynamics 
of  the  plant  it  observes  (Ref  6:190).  An  obvious  solution 
is  to  place  the  observer  poles  far  into  the  left  half  of 
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the  s plane,  but  this  introduces  differentiator  noise 
as  the  poles  become  very  large  negatively.  This  question 
will  be  investigated  with  the  goal  of  establishing  some 
guideline  for  selecting  observer  poles  to  minimize  system 
sensitivity  to  parameter  variations. 

The  question  of  what  constitutes  a sensitivity 
function  has  no  standard  answer.  Anderson  (Ref  10:117-124) 
considers  the  return  difference  matrix,  I+G (jw) H ( jw) , 
where  G(jw)H(jw)  is  the  open-loop  transfer  function,  as 
the  determining  factor.  Since  the  closed-loop  transfer 
function  is  G ( jw) / [I+G ( jw) H ( jw) ] , a system  obviously  is 
insensitive  to  a parameter  variation  in  the  plant,  G(jw), 
if  the  return  difference  is  "large."  On  the  other  hand, 
D'Azzo  and  Houpis  define  a sensitivity  function  to  be  the 
ratio  of  derivatives  of  the  natural  logarithms  of  the  con- 
trol ratio  and  the  variable  parameter.  Their  algebraically- 
reduced  sensitivity  function  is 


(11) 


where  S is  the  sensitivity  of  the  control  ratio,  M,  with 
respect  to  parameter  p.  Equation  (11)  is  a normalized  form 
of  the  more  general  sensitivity  function,  3y/3p,  (Ref  13:1) 
which  simply  expresses  the  change  in  output  due  to  a 
parameter  variation.  This  definition  follows  from  Eq  (10) 
and  will  be  used  as  the  standard  definition  of  sensitivity 
throughout  this  report. 
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Figure  3.  State-variable  Feedback 


Figure  3 depicts  the  application  of  state-variable 
feedback  to  a linear  system  that  can  be  modelled  by  a 
first-order  differential  equation.  The  Laplace  transform 
of  the  describing  differential  equation  for  the  system  is 

sX(s)  = AX (s)  + BU  (s)  (12) 

with  output  Y(s)  = CX(s)  (13) 

and  input  U(s)  * -KX(s)  + R(s)  (14) 

I 

Substituting  Eq  (14)  into  Eq  (12)  and  rearranging  terms 
yields 


sX(s)  = (A-BK)X(s)  + BR  (s)  (15) 


Taking  the  partial  derivative  of  Eqs  (15)  and  (13)  with 
respect  to  parameter  p yields,  respectively. 


3sX(s) 

“lp“ 


X(s)  + 


(A-BK) 


3X  (s) 

"5p"  + 


If 


(16) 


(17) 


L 


3Y(s) 

3p“ 


3C 


3X  (s) 


^X(s)  + C — 


It  is  noted  that  the  input  R(s)  is  not  a function  of  the 
variable  parameter.  The  partial  derivatives  are  evaluated 
at  the  nominal  parameter  value.  Eqs  (15)  and  (16)  may  be 
combined  to  form  an  augmented  set  of  matrix  equations,  or 


sX(s) 
3sX (s) 


3p 


A-BK 


L 


3 (A-BK) 

3p 


A-BK 


X(s) 
3X  (s) 


dp 


B 


3B 

3p. 


R(s) 


(18) 


3Y(s) 


3P 


! C 


3p 


X (s) 


3X  (s) 


L 3p  J 


(19) 


These  two  equations  have  the  same  form  as  Eqs  (12)  and 
(13) . For  simplicity,  the  2n  x 2n  augmented  "A"  matrix 
of  Eq  (18)  is  labelled  A'.  Then,  solving  Eq  (18)  for  the 
augmented  "X  (s) " vector  and  substituting  it  into  Eq  (19) 
yields  the  sensitivity  transfer  function 


3Y(s) 

5p“ 


3C  i „ 
rr1  ' C 
3p  l 

l 


sI-A' 


-1 


B 


3B 

3P. 


R(s) 


(20) 
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In  Fig.  4 an  observer  has  been  added  to  the  basic 
system  of  Fig.  3 to  provide  the  reconstructed  state  vector, 
X(s) . An  acceptable  choice  for  T is  the  identity  matrix 
so  that  the  state  vector  of  the  observer  is,  in  steady 
state,  the  state  vector  of  the  plant  it  is  observing. 


Figure  4 . Use  of  Observer  in  State-variable  Feedback 


Choosing  T=I  also  reduces  the  mathematics  involved  since 
F»B  and  Eq  (8)  becomes 

A-D=VC  (21) 


In  the  frequency  domain,  with  the  substitution  U(s)= 
-K(s)+R(s)  , Eq  (7)  can  be  expressed  as 

sZ(s)  - (D-BK)Z(s)  + VCX(S)  + BR  (s)  (22) 

Considering  that  the  matrices  D and  V are  not  functions 
of  the  plant  parameter  p,  since  they  are  computed  at  the 
nominal  value  of  p,  the  partial  derivative  of  Eq  (22)  with 


respect  to  p is 


3sZ (s) 


KZ(s)  + (D-BK) 


3Z  (s) 


+ v § X(s) 


3X(s) 


+ VC 


+ If  *(s) 


As  in  the  state-variable  feedback  case  it  is  desired  to 
construct  an  augmented  matrix  equation  as  a prelude  to  the 
sensitivity  transfer  function  for  the  system  in  Fig.  4. 
Combining  Eas  (12),  (22),  (16),  and  (23)  gives 


sX  (s) 
sZ^(s) 

3sX (s) 


3sZ  (s) 
" 9P~ 


-BK 

D-BK 


V— 

3p 


X (s) 
Z(s) 

I 3 X (s) 


D-BK  - 


3Z  (s) 


+ 3B  R(s) 
3P 

3B 

3p 


As  before,  Eq  (24)  can  be  solved  for  the  augmented  state 
vector,  which  is  substituted  into  Eq  (17)  to  obtain  the 
sensitivity  transfer  function.  Defining  the  partitioned 
4n  x 4n  "A"  matrix  of  Eq  (24)  as  A",  the  final  result  is 


3Y(s) 


B 

-1  __B__ 


0 ! C ol  fsI-A,'l  3B  R(s) 

L3P  ! JL  J TB  ~ 


where  each  zero  element  in  the  row  matrix  above  is  itself 
an  rxn  matrix  of  zeros  to  maintain  conformability . 
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This  concludes  the  theoretical  development  required 
for  the  problem.  Since  it  is  obviously  not  practical  to 
analytically  compare  Eq  (20)  with  Eq  (25)  to  determine  if 
one  system  is  more  sensitive  to  a parameter  variation  than 
the  other,  the  next  chapter  will  develop  numerous  examples 
to  investigate  this  question. 


III.  Numerical  Examples 

In  order  to  keep  the  required  computations  to  a 
convenient  level  a two-state  plant  with  one  input  and  out- 
put will  be  examined.  Even  such  a simple  example  will 
give  some  insight  into  the  two  sensitivity  transfer  func- 
tions and,  more  importantly,  an  idea  of  how  the  three 
areas  under  investigation  affect  the  sensitivity  of  a 
system. 

Consider  the  plant 

in  Fig.  5 with  nominal  param-  

u _k_  i 

eter  values  of  K=10,  a=2,  and  | s+a  | s+b 

b“l  where  "a"  is  the  varia-  

ble  parameter.  The  states  Figure  5.  Example  Plant 

i Xi  and  X2  will  be  defined  as  outputs  of  the  integrators 
from  right  to  left  or  top  to  bottom,  respectively,  as 
applicable  for  each  situation. 

In  physical  variables  the  state  equations,  omitting 
the  time  designators  for  brevity,  are 


Figure  5.  Example  Plant 


-11  0 

x + u 
0 -a  ~ 10 


[l  o]  X 


The  implementation  of  state-variable  feedback  is  shown  in 
Laplace- transformed  form  in  Fig.  6.  Using  block  diagram 
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Figure  6.  Physical  Variable  Feedback 
reduction  techniques  the  control  ratio  is  found  to  be 


a 

r s2+(3+10k2)s+(2+10k2+10ki) 

To  generate  some  sort  of  performance  specifications,  sup- 
pose it  is  desired  that  the  system  have  zero  steady-state 
error  for  a step  input.  Then  2+10k2+10k j=10  and  the  con- 
trol ratio  is 

1=  12 (30) 

s!+(3+10ki)s+10 

Now  the  pole  placement  is  a function  only  of  k2  with  kj 
being  determined  by  2+10k2+10ki=10.  Four  cases  of  pole 
placement  will  be  investigated.  Hereafter,  "plant  poles" 
are  considered  to  be  the  closed-loop  poles  of  the  plant. 


"" 


Case  1.  Plant  poles:  -2,  -5;  physical  variables. 

The  plant  poles  will  be  shifted  from  -1,  -2  to 
-2,  -5  to  provide  faster  transient  response  decay.  The 
characteristic  equation  is 

(s+2)  (s+5)  = s2  + 7s  + 10  = 0 (31) 

Equating  the  plant  characteristic  equation  (denominator 
of  Eq  (30))  with  Eq  (31)  yields  k2=0.4  and  thus,  k1=0.4. 

Case  2.  Plant  poles:  -2 . 2 35+ j2 . 235;  physical  variables. 

The  plant  poles  will  be  placed  at  -2 . 235+ j 2 . 235 
to  give  a slightly  oscillatory  response  with  a damping 
ratio  0f  0.707.  The  characteristic  equation  is 

(s+2. 235+ j 2. 235)  (s+2 . 2 35- j 2 . 235 ) = s2  + 4 . 4 7s  + 10  = 0 (32) 

with  k2  = 0.1472  and  k2  = 0.6528  . 

Case  3.  Plant  poles:  -1,  -10;  physical  variables. 

For  a larger  shift  to  the  left,  poles  are  placed 
at  -1,  -10.  Because  of  its  rapid  transient  decay,  the 
pole  at  -10  is  considered  to  be  nondominant.  The  charac- 
teristic equation  becomes 

(s+1)  (s+10 ) = s2  + lls  + 10  = 0 (33) 

making  k2  = 0.8  and  kj  = 0.0. 

Case  4.  Plant  poles:  -3.125,  -3.2;  physical  variables. 

Finally,  the  two  poles  are  placed  close  together 
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at  -3.125  and  -3.2.  According  to  Singer  (Rof  4:88),  this 
situation  has  a high  eigenvalue  sensitivity.  This  test 
will  help  determine  the  extent,  if  any,  eigenvalue 
sonsitivity  affects  system  sensitivity.  The  resulting 
characteristic  equation  is 

(s+3.125)  (s+3.2-  s2+  6.325s  + 10  = 0 (34) 

forcing  k2  = 0.3325  and  k)  = 0.4675. 

Three  observer  designs  will  bo  considered,  each 
having  judiciously-placed  poles.  It  is  hoped  that  this 
exercise  will  provide  some  guidance  for  observer  poio 
selection  for  minimizing  sensitivity  to  plant  parameter 
variations. 


Observer  I.  Poles:  -6,  -6;  physical  variables. 

Double  polos  are  arbitrarily  selected  at  -6.  For 
all  but  Case  3 above  this  choice  places  the  observer  poles 


to  the  left  of  the  plant  poles  so  that  the  observer 
reaches  steady  state  before  the  plant.  Analysis  of  Case  3 
results  should  suggest  whether  or  not  this  requirement  is 
valid.  With  a characteristic  polynomial  of  s'+12s+36  and 
T-I,  the  observer  dynamics  matrix,  from  Eq  (21),  is 


"3  1 1 d : 2 -1-v  i 1 

3 2 J d j 2 “V2  -2 


(35) 


The  characteristic  polynomial  of  D is  found  by  equating  the 
determinant  of  [gl-D]  to  zero,  or 


1 
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s?  + (3  + \'i  ) s+  (2+2v,  + Vj ) - 0 


(36) 


Equating  the  coefficients  of  s in  Eq  (36)  to  those  of  the 
desired  characteristic  polynomial,  s7+12s+36,  V and  D are 
found  to  be 


9 

16 


(37) 


D, 


-10  1 
-16  -2 


(38) 


Observer  II.  Poles:  -12,  -15;  physical  variables. 

By  selecting  poles  to  the  left  of  -10,  the  observer 
dynamics  will  always  have  reached  steady  state  before  the 
plant.  A choice  of  -12,  -15  accomplishes  this  objective 
plus  it  provides  a large  separation  between  plant  and 
observer  poles  for  Cases  1,  2,  and  4.  With  a characteris- 
tic polynomial  of  s'+27s+180,  V and  P are  computed,  as 
before  with  Observer  I,  to  be 


V,  = 


24 

130 


(39) 


-25  1 

-130  -2 


(40) 


Observer  III.  Poles:  -3,  -4;  physical  variables. 

For  the  third  trial,  observer  poles  will  be 
located  in  tho  vicinity  of  the  plant  poles  to  allow  some 
interaction  between  the  two  transient  responses.  Roots 
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A 


of  -3,  -4  generate  a characteristic  polynomial  of 
s2+7s+12,  making 


V 

s 


(41) 


D ** 

3 


1 

-2 


(42) 


In  phase-variable  form  the  state  equations  for  the 
plant  in  Fig.  5 are  derived  from  the  describing  differen- 
tial equation  to  be 


0 

1 

x + 

r°i 

-a 

-a-1 

L i J 

(43) 


y - [10  0]  x (44) 

The  block  diagram  for  this 
system  is  given  in  Fig.  7. 

The  variables  to  be  fed 
back  are  not  the  same  xj 
and  x 2 as  before;  there-  Figure  7.  Phase  Variable  Plant 
fore,  a new  relationship  b 'tween  the  control  ratio  and 
feedback  coefficients  is  necessary.  Equation  (15)  can  be 
solved  for  X which  is  substituted  into  Fq  (13)  to  yield 
the  general  matrix  form  of  the  control  ratio 


£ - C[sI-A+BK]“‘b 


(45) 


I 


Substituting  the  corresponding  matrices  of  Eqs  (431  and 
(44)  into  Eg  (45)  yields 


r 


10 

b2+  (k2+3)  s+  0^+2) 


(46) 


which  is  the  desired  control  ratio  as  a function  of  phase- 
variable  feedback  coefficients. 

The  necessary  feedback  coefficients  are  to  be 
calculated  for  pole  placements  identical  to  those  presented 
in  the  physical  variable  development  using  the  same  pro- 
cedure discussed  there.  It  is  obvious  that  for  zero 
steady-state  step  error  ki=8  and  is  independent  of  pole 
placement. 

Case  5.  ^lant  poles:  -2,  -5;  phase  variables. 

The  characteristic  equation  was  s2+7s+10=0,  so 
that  k2  *4 . 

Case  6.  Plant  poles:  -2.235+j  2.235;  phase  variables.  ! 

I 

The  characteristic  equation  was  s2  + 4 . 471s+10=0 , 
making  k2 *1.471. 

Case  7.  Plant  poles:  -1 , -10;  phase  variables. 

The  characteristic  equation  was  s2+lls+10=0,  or 

j 

k2-8. 

Case  8.  Plant  poles:  -3.125,  -3.2;  phase  variables. 

The  characteristic  equation  was  s2+6 .325s+101“0 , 

thus  k2"3.325. 

pi 

' 

22 

■ 

1 1 


i 
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Observer  IV.  Poles:  -6,  -6;  phase  variables. 

Solving  Eq  (21)  at  nominal  parameter  values  and 
equating  the  coefficients  of  the  characteristic  poly- 
nomial of  D with  the  coefficients  of  the  desired  character- 
istic equation,  s*+12s+36=0,  results  in 


Vh 


[:::] 

E;  -;l 


(47) 


(48) 


Observer  V.  Poles:  -12,  -15;  phase  variables. 

Repeating  the  procedures  for  determining  V and  D 
with  the  desired  observer  characteristic  equation  of 
s*+27s+180®0  gives 


M 

[10.6J 


[- 


-24  1 

108  -3 


(49) 


(50) 


Observer  VI.  Poles:  -3,  -4;  phase  variables. 

For  the  desired  characteristic  equation  of 
sa+7s+12*0,  V and  D were  computed,  as  before,  to  be 


[■: .;] 


(51) 


(52) 
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To  describe  the  basic  plant  in  Jordan  canonical  form,  the 
plant  transfer  function  must  be  rewritten  in  its  partial 
fraction  expansion  form  as 


W (_1 1 ) u 

1-a  's+a  s+l; 


The  simulation  diagram 
for  Eq  (53)  is  given 
in  Fig.  8.  The  new 
state  variables  X]  and 
X2  are  not  the  same  as 
in  either  of  the  pre- 
ceding state-variable 
forms.  The  state 
equations  can  be 
derived  from  Fig.  8 


Figure  8.  Jordan  Form  Plant 


to  be 


t: .:]  **[:]" 

! 10  -10  v 
IT^T  T^aJ  - 


To  obtain  the  control  ratio  as  a function  of  the  feedback 
coefficients  in  this  case,  the  A,  B,  and  C matrices  of 
Eqs  (54)  and  (55)  are  substituted  into  Eq  (45) . The  result 


r sa+(3+k1+k2)s+(2+kI+2k2) 


For  zero  steady-state  step  error,  2+kj+2k2=10.  The  feed- 
back coefficients  are  computed  for  each  case  by  equating 
like  coefficients  of  s in  Eq  (56)  with  those  of  the 
desired  characteristic  equations. 

Case  9.  Plant  poles:  -2,  -5;  Jordan  canonical  variables. 

The  characteristic  equation  was  sz+7s+10=0; 
ki-0,  k2=4 . 

Case  10.  Plant  poles:  -2.235+j2 .235;  Jordan  canonical 
variables. 

The  characteristic  equation  was  s2  + 4 .471s+10=0; 
kj— 5.056,  kz  = 6.528. 

Case  11.  Plant  poles:  -1,  -10;  Jordan  canonical  variables. 
The  characteristic  equation  was  s2+lls+10=0; 

ki*8 , k 2 = 0 • 

Case  12.  Plant  poles:  -3.125,  -3.2;  Jordan  canonical 
variables. 

The  characteristic  equation  was  s2+6.325s+10=0; 
ki*-l .35,  k2=4 . 675. 

Using  the  same  procedure  described  in  the  physical 
and  phase  variable  cases,  the  following  V and  D matrices 
were  computed  for  the  Jordan  form  observers: 


Observer  VII.  Poles:  -6,  -6;  Jordan  canonical  variables 


T7 T 


” 7 — e ' 


mm 


.. 


4, 


C57) 


(58) 


Observer  VIII.  Poles:  -12,  -15;  Jordan  canonical  variables. 


(59) 


(60) 


Observer  IX.  Poles:  -3,  -4;  Jordan  canonical  variables. 


(61) 


(62) 


Sensitivity  transfer  functions  were  computed  for 
each  of  the  48  combinations  of  state-variable  representa- 
tions, pole  placements,  and  observer  designs  as  displayed 
in  Fig.  9.  The  algorithm  used  to  perform  the  necessary 
matrix  inversions  for  Eqs  (20)  and  (25)  is  presented  in 
Appendix  A.  As  a means  for  comparing  relative  sensitivi- 
ties of  the  state-variable  feedback  model  with  the 
observer  feedback  models,  frequency  responses  of  the 
sensitivity  transfer  functions  were  plotted  from  0-5 
radians  per  second.  This  range  is  considered  to  be 
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representative  of  the  operational  bandwidth  of  practical 
control  systems.  Figs.  10-17  display  the  absolute  value 
of  the  magnitudes  of  the  sensitivity  functions  for  all 
the  cases.  It  was  noticed  that  the  results  for  physical 
variables  and  phase  variables  were  always  identical  (to 
about  four  decimal  places);  therefore,  only  one  plot  was 
used  to  represent  the  two  data  runs  for  each  similar  test. 


Plant  Physical  Phase  Jordan  Canonical 

Poles  Variables  Variables  Variables 


'Tl! 
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Figure  9.  Matrix  of  Test  Cases 


SENSITIVITY  MAGNITUDE 


FREQUENCY -RADIANS  PER  SECOND 

0 Direct  state  feedback. 

O Observers  I and  IV.  Poles:  -6,-6. 

□ Observers  II  and  V.  Poles:  -12,-15. 

A Observers  III  and  VI.  Poles:  -3,-4. 


Figure  10.  Case  1.  Plant  Poles:  -2,-5;  Physical  Variables. 

Case  5.  Plant  Poles:  -2,-5;  Phase  Variables. 
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SENSITIVITY  MAGNITUDE 


! 


FREQUENCY -RADIANS  PER  SECOND 


0 Direct  state  feedback. 

O Observer  VII.  Poles:  -6,-6. 

□ Observer  VIII.  Poles:  -12,-15. 
A Observer  IX.  Poles:  -3,-4. 


Figure  11.  Case  9.  Plant  Poles:  -2,-5;  Jordan  Canonical  Variables. 


30 


SENSITIVITY  MAGNITUDE 


FREQUENCY-RADIANS  FER  SECOND 

0 Direct  state  feedback. 

O Observer  VII.  Poles:  -6,-6. 

□ Observer  VIII.  Poles:  -12,-15. 

A Observer  IX.  Poles:  -3,-4. 


Figure  13.  Case  10.  Plant  Poles:  -2.235±j2.235;  Jordan  Canonical 
Variables . 


FREQUENCY -RADIANS  PER  SECOND 

0 Direct  state  feedback. 

O Observers  I and  IV.  Poles:  -6,-b. 

□ Observers  II  and  V.  Poles:  -12,-15. 

A Observers  III  and  VI.  Poles:  -3»-4. 


Figure  14.  Case  3.  Plant  Poles:  -1,-10;  Physical  Variables 
Case  7.  Plant  Polos:  -1,-10}  Phase  Variables. 


SENSITIVITY  MAGNITUDE 


0 Direct  state  feedback. 

O Observer  VII.  Poles:  -6,-6. 

□ Observer  VIII.  Poles:  -12,-15. 
A Observer  IX.  Poles:  -3,-4. 


Figure  15.  Case  11.  Plant  Poles:  -1,-10}  Jordan  Canonical  Variables. 
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SENSITIVITY  MAGNITUDE 


6 

3 


1 

.6 

.3 


.1 

0 12  3 4 5 

FREQUENCY -RADIANS  FER  SECOND 

Direct  state  feedback. 

O Observers  I and  IV.  Poles:  -6,-6. 

□ Observers  II  and  V.  Poles:  -12,-15. 

A Observers  III  and  VI.  Poles:  -3,-4. 

Figure  16.  Case  4«  Plant  Poles:  -3.125,-3.2;  Physical  Variables. 

Case  8.  Plant  Poles:  -3.125,-3.2;  Phase  Variables. 


FREQUENCY -RADI ANS  FER  SECOND 

0 Direct  state  feedback. 

O -Observer  VII.  Poles:  -6,-6. 

□ Observer  VIII.  Poles:  -12,-15. 

A Observer  IX.  Poles:  -3,-4* 


Figure  17.  Case  12.  Plant  Poles:  -3.125,-3.2;  Jordan  Canonical 
Variables . 


IV.  Results  and  Discussion 


To  avoid  confusion  over  notation  of  the  different 
test  cases  in  this  section,  an  explicit  coding  is  needed  to 
distinguish  the  three  aspects  of  the  cases:  state-variable 
representation,  plant  poles,  and  observer  pole  locations. 
Since  the  results  for  physical  variables  and  phase  vari- 
ables were  identical  in  all  cases  studied,  it  is  necessary 
to  refer  only  to  one  of  the  two;  say,  phase  variables.  Then 
the  following  notation  is  adopted:  (a;  b,  c;  d,  e)  where  a 
is  the  type  of  state  variable  (P  for  phase  variables?  J for 
'Jordan  canonical  variables) ; b and  c are  the  plant  pole 
locations  (for  complex  conjugate  poles  b and  c are  the  real 
and  imaginary  parts,  respectively) ; d and  e are  the  observer 
pole  locations.  If  d and  e are  not  listed,  direct  state 
feedback  is  implied.  When  observers  are  individually  men- 
tioned by  their  respective  numbers,  the  type  of  state  vari- 
able representation  and  observer  poles  will  immediately 
follow  as  in  IV  (P;  -6,  -6). 


For  the  phase-variable  representation  (Fig.  10),  all  three 
observers  exhibited  greater  sensitivity  magnitudes  than  did 
the  plant  with  direct  state  feedback;  in  Jordan  form 
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(Fig.  11)  the  opposite  was  true  with  an  increase  of  two 
orders  of  magnitude  in  the  direct  state  feedback.  The 
observer  sensitivities  were  independent  of  state-variable 
representation.  Although  observers  V (P;  -12,  -15)  and 
VIII  (J;  -12,  -15)  were  the  least  sensitive,  followed  by 
observers  VI  (P;  -3,  -4)  and  IX  (J;  -3,  -4),  then  IV 
(P;  -6,  -6)  and  VII  ( J ; -6,  -6)  in  increasing  order,  none 
had  a sensitivity  magnitude  greater  than  0.54  in  the  band- 
width of  interest;  thus  the  system  changes  by  only  half  the 
magnitude  of  a parameter  variation.  Observers  IV  (P;  -6, 

-6),  VI  (P?  -3,  -4),  VII  (J;  -6,  -6)  and  IX  (J?  -3,  -4) 
displayed  extreme  insensitivity  to  parameter  variations  at 
about  one  radian  per  second;  observers  V (P;  -12,  -15)  and 
VIII  (J;  -12,  -15)  were  insensitive  near  two  radians  per 
second,  as  was  the  direct  state  feedback  function  in  phase- 
variable  form.  In  Jordan  form  the  direct  state  feedback 
sensitivity  was  not  affected  by  this  phenomenon.  The 
algorithm  used  to  generate  data  points  considered  only  the 
absolute  value  of  the  sensitivity  function.  In  expanding 
the  [sI-A] ~ 1 terms,  some  of  the  components  were  suffi- 
ciently negative  to  impart  a negative  sign  to  the  sensitivity 
function  for  some  frequencies.  The  physical  implication  is 
that  the  system  output  deviates  in  the  opposite  direction 
of  the  parameter  variation.  It  also  means  there  is  a fre- 
quency for  which  the  system  has  zero  sensitivity  to  param- 
eter variations. 
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Case  6.  (P?  -2.235,  ±2.235),  (P;  -2.23S,  ±2.235; 
-6,  -6),  (P;  -2.235,  +2.235?  -12,  -15),  (P;  -2.235,  +2.235; 


-3,  -4)  and  Case  10.  (J;  -2.235,  +2.235),  (J;  -2.235, 
+2.235?  -6,  -6),  CJ;  -2.235  +2.235;  -12,  -15),  (J?  -2.235, 


+2.235;  -3,  -4) . In  Fig.  12  observers  IV  CP;  “6,  -6)  and 
V (P;  -12,  -15)  very  closely  approximated  the  sensitivity 
function  for  direct  feedback.  Observer  VI  (P?  -3,  -4) 
exhibited  the  sign  inversion,  previously  discussed,  at  about 
0.3  radians  per  second,  becoming  extremely  sensitive  with 
increasing  frequency.  In  Fig.  13,  observer  IX  CJ;  -3,  -4) 
was  even  more  sensitive  to  variations,  being  over  two 
orders  of  magnitude  above  observers  VII  (J;  -6,  -6)  and 
VIII  (J;  -12,  -15).  Again,  in  Jordan  form,  the  direct  state 
feedback  system  sensitivity  was  an  order  of  magnitude 
greater  than  it  was  in  phase-variable  form.  Observers  IV 
(P;  -6,  -6),  V (P?  -12,  -15),  VII  (J;  -6,  -6),  and  VIII 
(J;  -12,  -15)  were  all  very  similar  with  peak  sensitivities 
of  0.3. 

Case  7.  (P;  -1,  -10),  (P;  -1,  -10;  -6,  -6), 

-3,  -4)  and  Case  11. 


-12,  -15) 


(J;  -1,  -10;  -3,  -4) . The  only  difference  between  Figs.  14 
and  15  was  the  sensitivity  of  the  direct  state  feedback 
system.  In  phase-variable  form,  direct  state  feedback  was 
the  least  sensitive  implementation  whereas  in  Jordan  form, 
it  was  increased  by  two  orders  of  magnitude;  observer 
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performances  were  identical  for  the  two  state-variable 
representations . 

Case  8.  (P;  -3.125,  -3.2),  (P;  -3.125,  -3.2; 

-6,  -6),  (P;  -3.125/  -3.2;  -12,  -15),  (P;  -3.125,  -3.2; 

-3,  -4)  and  Case  12.  (J;  -3.125,  -3.2),  (J;  -3.125,  -3.2; 

-6,  -6),  (J;  -3.125,  -3.2;  -12,  -15),  (J;  -3.125,  -3.2; 

-3 , -4) . In  these  cases  all  systems  exhibited  very  high 
sensitivities  to  parameter  variations  with  the  observers 
being  basically  independent  of  state-variable  representa- 
tion. Observer  VI  (P;  -3,  -4)  underwent  a sign  inversion 
at  0.2  radian  per  second;  observers  VII  (J;  -6,  -6)  and 
VIII  (J;  -12,  -15)  experienced  sign  inversions  at  0.55  and 
0.8  radian  per  second,  respectively.  The  direct  state  feed- 
back system  sensitivity  was  much  greater  in  Jordan  form 
than  in  phase-variable  form. 

Analysis 

The  first  question  posed  in  developing  this  study 
concerned  the  effects  of  different  state-variable  representa- 
tions upon  system  sensitivity.  The  data  show  that  there  is 
no  difference  in  sensitivity  between  physical  and  phase 
variables  and  that  the  Jordan  form  is  always  more  sensitive 
than  the  former  two  representations.  This  does  not  corre- 
late with  eigenvalue  sensitivity  concepts.  Reference  4 
demonstrates  that  unity  eigenvalue  sensitivity  is  always 
obtained  with  the  Jordan  form  if  the  eigenvalues  are  real 
whereas,  with  phase  variables,  the  eigenvalue  sensitivity  can 
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range  from  less  than  one  to  infinity,  depending  upon  the 
pole  locations.  The  following  phase-variable  eigenvalue 
sensitivities  (S)  were  computed  for  the  examples  investi- 
gated: for  plant  poles  of  -2,  -5:  S=l,  2,  respectively; 
for  plant  poles  of  -2.235  + j 2.235:  S=0.88,  0.88, 
respectively;  for  plant  poles  of  -1,  -10:  S=0.22,  1.22, 
respectively;  and  finally,  for  plant  poles  of  -3.125,  -3.2: 
S»55,  56,  respectively  (Ref  4:88).  The  phase-variable 
eigenvalue  sensitivities  for  the  first  three  sets  of  plant 
poles  are  all  in  the  vicinity  of  S=l,  which  is  the  Jordan 
form  eigenvalue  sensitivity,  and  yet  the  system  sensitivi- 
ties for  those  cases  vary  greatly  between  the  two  different 
state-variable  representations.  It  is  not  understood  why 
this  apparent  discrepancy  occurs.  It  would  seem  that  low 
system  sensitivity  would  imply  low  eigenvalue  sensitivity 
and  conversely,  but  the  converse  does  not  hold  here.  The 
observer  sensitivities,  on  the  other  hand,  did  not  vary  much 
from  one  state-variable  representation  to  another;  only  the 
direct  state  feedback  sensitivity  functions  were  affected 
by  the  Jordan  form. 

The  second  question  was  that  of  pole  placement 
effects  upon  system  sensitivity.  The  extent  of  desired 
pole  placement  definitely  can  change  the  system  sensitivity. 
Attempting  to  place  poles  close  together  drastically 
increases  parameter  variation  sensitivity;  creation  of 
complex  conjugate  poles  does  not  have  adverse  sensitivity 
effects.  Low  system  sensitivity  is  also  achieved  by  pairing 
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a pole  near  the  origin  with  one  deeper  into  the  left  half- 
plane. The  magnitudes  of  the  feedback  gains  required  to 
achieve  the  pole  placements  do  not  appear  to  affect  the 
system  sensitivities  as  seen  by  the  duplicate  results  for 
physical  and  phase  variables  for  which  the  feedback 
coefficients  differed  by  an  order  of  magnitude.  There 
were  no  obvious  adverse  effects  from  the  negative  ki 
coefficients  (which  create  positive,  unstable  feedback) 
in  several  of  the  Jordan  form  cases  as  all  exhibited  the 
same  high  sensitivity  characteristics. 

The  final  area  of  investigation,  that  of  effects  of 
observer  design  upon  system  sensitivity,  revealed  a definite 
relationship  between  the  two  factors.  The  design  of 
observer  dynamics  is  the  most  important  aspect  in  reducing 
the  sensitivity  of  an  observed  plant  to  internal  parameter 
variations.  As  mentioned  in  Chapter  II,  literature  on 
observers  states  only  that  the  observer  poles  should  be 
placed  to  the  left  of  all  plant  poles  for  good  performance. 
Analysis  of  the  test  cases  in  this  report  indicates  that 
such  a guideline  may  be  too  restrictive  upon  observer  design. 
In  several  of  the  cases  involving  plant  poles  of  -1,  -10, 
placing  observer  poles  to  the  left  of  the  dominant  pole  at 
-1 , but  to  the  right  of  the  nondominant  pole  at  -10,  achieved 
not  only  low  system  sensitivity  but  identical  performance  in 
all  three  state-variable  forms.  In  general,  though,  placing 
the  observer  poles  to  the  left  of  all  plant  poles  does 
result  in  a lower  system  sensitivity. 
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V.  Conclusions  and  Recommendations 


Conclusions 

The  selection  of  state  variables  to  be  utilized  in 
system  analysis  does  affect  the  outcome.  Physical  vari- 
ables are  the  best  choice  when  applicable  because  they 
represent  measurable  quantities  of  interest  within  a system. 

If  physical  variables  are  not  realizable  for  some  reason, 
such  as  when  a derivative  of  the  input  (zero)  appears  with 
the  n-th  state,  phase  variables  yield  identical  sensitivity 
results.  Jordan  canonical  variables  should  be  avoided 
because  of  the  relatively  high  system  sensitivity  that 
results  from  this  form. 

In  placing  closed-loop  poles  through  state-variable 
feedback,  system  sensitivity  is  reduced  when  the  poles  are 
separated  from  each  other,  especially  if  one  pole  is  non- 
dominant. This  result  is  also  seen  in  Chapter  12  of  Reference 
1.  Observer  poles  should  be  placed  to  the  left  of  all  dominant 
plant  poles  as  a minimal  requirement,  and  to  the  left  of  all 
plant  poles,  if  practical,  for  low  system  sensitivity  to 
parameter  variations. 

The  "bottom  line"  is  that  the  observer  can  be 
designed  so  that  the  resulting  system  sensitivity  is  as  low 
as  that  of  direct  state  feedback,  or  even  lower  in  many  cases. 

Recommendations 

Several  areas  for  further  research  were  identified 
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during  this  study.  The  observers  considered  were  full- 
order  types;  that  is,  the  entire  state  vector  was  recon- 
structed for  feedback  purposes.  It  is  possible  to  use 
available  states  for  feedback  and  design  a reduced-order 
observer  to  generate  only  the  inaccessible  states.  The 
observer  has  fewer  states  and  may  reduce  system  sensitivity 
since  there  is  less  probability  of  interaction  between  plant 
parameter  variations  and  the  observer.  There  exists  a 
trade-off  between  acceptable  system  sensitivity  and  the 
cost/availability  of  sensors  to  pick  off  plant  states.  An 
investigation  should  be  conducted  to  determine  the  sensi- 
tivity aspects  of  reduced-order  observers. 

Equation  (10)  expresses  the  change  in  output  magni- 
tude due  to  a parameter  variation,  but  does  not  address 
corresponding  phase  angle  changes.  In  the  data  summary  of 
Chapter  IV,  sign  changes  were  credited  with  causing  zero 
parameter  sensitivities  at  certain  frequencies.  A question 
arises  as  to  whether  or  not  it  is  possible  to  predict  condi- 
tions under  which  the  "phase  angle  sensitivity"  (3<#>/3p) 
causes  a sign  change  in  the  sensitivity  magnitude.  If  so, 
it  may  be  possible  to  alter  some  design  parameters  to  create 
zero  sensitivity  at  frequencies  of  known  external  distur- 
bances. Such  a capability  would  be  a most  useful  addition 
to  control  theory. 

Another  area  is  that  of  plants  with  multiple  inputs 
and  outputs.  Luenberger  states  that  the  problem  of  design- 
ing an  observer  for  a plant  with  multiple  outputs  involves 


a series  of  observers,  one  for  each  output  (Bef  6:190). 
Since  the  output  of  each  observer  is  intimately  tied  to  the 
plant,  and  hence  to  each  other  observer  through  cross- 
coupling, there  may  be  some  interesting  sensitivity  conse- 
quences to  a plant  parameter  variation. 
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The  only  difficult  task  in  computing  the  sensi- 
tivity transfer  functions,  Eqs  (20)  and  (25),  is  inverting 
the  2n  x 2n  matrix  [sI-A']  and  the  4n  x 4n  matrix  [sI-A"J • 
Obviously,  for  n>l  this  becomes  a formidable  task  for  hand 
calculation  and  computer  routines  inherently  solve  numeri- 
cal matrix  inversions.  Reference  11  presents  a convenient 
scheme  for  accomplishing  this  task,  based  upon  the  fact 
that 

, [sI-M]  “ 1 = *(eMt)  (63) 


There  are  a number  of  methods  available  for  expanding  a 
matrix  exponential,  but  this  scheme  is  readily  adaptable 
to  a digital  computer.  The  four  basic  steps,  each  of 
which  is  discussed  later,  are  as  follows: 

1.  Construct  the  Vandermonde  matrix,  A. 

2.  Compute  A-1  - V. 

3.  Construct  the  component  matrices,  of  M. 

4.  Expand  Eq  (63)  as 


X (eMt) 


q "hr1 
- Z E 
k-1  i-0 


Jk,i 


i ^k^ 
*(tAe  k ) 


(64) 
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where  q is  the  number  of  distinct  eigenvalues  of  M and 


"k is  the  multiplicity  of  in  the  minimal  polynomial  of 
M.  (See  Appendix  B for  a brief  discussion  of  minimal 
polynomials. ) 

A useful  property  of  the  matrices  A'  and  A"  in 
Eqs  (20)  and  (25) , respectively,  is  that  they  have  the 


form 


- fM  • °" 

M * H 


Jp  nj2n  x 2n 


where  n is  the  dimension  of  each  partition  of  M.  The 
eigenvalues  of  M are  simply  the  eigenvalues  of  M,  each 
repeated  once.  The  characteristic  polynomial  of  M is 
defined  as 


A (A)  - det  (AI-M)  = n ( A— A.  ) 2nk 

k-1  K 


The  Vandermonde  matrix  of  M is  constructed  from  the 
eigenvalues  of  M as  r rows  of  the  form 


(1  Ak  Ak*  ...  Ak2n"1] 


where  k»l,2,...,r  and  j=0,l. . . ,2nk~l  for  r real,  distinct 
eigenvalues  and  q-r  rows  of  the  form 


Real  (-^-j  [1  Ak  Ak*...  A^11"1] ) 


Iroag  (-^  -j  [1  Ak  X^2. . . X^2*1  XJ) 


where  k«r+l,  r+3,  ...  q-1  and  j“0,l, ... , 2nk-l  for  complex 
eigenvalues.  The  odd  subscript  accounts  for  the  fact  that 
complex  eigenvalues  occur  in  conjugate  pairs  as 


Xk  " °k+^wk'  Xk+1  “ ak"^wk 


Once  the  Vandermonde  matrix  is  complete.  Step  2 can  be 
accomplished  by  any  matrix  inversion  routine. 

Computation  of  the  components,  Z.  . of  M 

K 9 i 9 

depends  intimately  upon  the  form  of  the  Vandermonde  matrix 
for  each  system.  ^ is,  itself,  a partitioned  matrix  of 

nxn  component  matrices,  i and  ZPk  . of  M and  3M/3p, 
respectively;  that  is, 


where  Z, 


and  ZP. 


2n  -ui 
E M3  1 

j-1 


k-0,l,...q;  i'O,!,...,!^-] 


.u  j“l 

& 'If'  vj.*p 


k-0,1, 


rq;  i"0 , 1 , . . .n.-l 


The  subscript  fp  in  Eq  (71)  acts  as  a counting  index.  For 
example,  in  ZP^  Q fp«l;  ZP^  ^ fp=2  and  so  forth  up  to 
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fp««4n.  It  should  be  noted  that  fp  remains  fixed  over  the 
entire  summation  for  each  component  matrix  ZP^  The 
subscript  f in  Eq  (70)  is  not  so  simply  stated.  It  is  the 
index  number  of  the  row  in  the  Vandermonde  matrix  that 
corresponds  to  the  i-th  derivative  of  the  k-th  eigenvalue. 
This  is  better  understood  through  an  example.  Assume 
n«4  and  that  there  are  four  distinct,  real  eigenvalues  of 
M.  Then  q«4,  r*4,  and  ni*  n2*  nj*  n^  1.  The  Vandermonde 
matrix  is 


A ■ 


Ai6  Aj7 

6 A,5  7A  i* 

01  2A2  ... 

1 A,  A,2  ... 

0 1 2Aj  ... 

1 A^  Ah2  ... 

01  2Ah  ... 


1 A i Ai2  Ai*  A i*  A i* 

0 1 2A,  3A  i2  4 Ax5  SA^ 


(72) 


In  this  case,  for  Z ^ Q-*-f*l  ^the  0_th  derivative  of  \x  is 
row  1 of  A);  Z2^0^f«3;  Z3  g-*-f»5;  Z4  Q-*-f»7.  If,  say,  Ai, 
is  repeated  once,  then  q»3=*r  and  ni«2,  nj-nj-l  and  the 
Vandermonde  matrix  is 
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F°r  z,  Q-^f=l;  Z,  ^->-f=2  (the  first  derivative  of  is  row 


2 of  A) ; Z2  Q-*-f=5;  Z^  Q-*-f=7.  The  pattern  holds  for  any 
case.  Once  Z,  . and  ZP.  . have  been  computed,  the  component 

JC  / 1 K f 1 

matrices  Z,  . are  formed  as  shown  in  Eq  (69) . It  should  be 
noted  that  for  each  Zv  . there  is  a ZPV  . but  not  neces- 
sarily  a Z^  ..  In  such  cases  the  missing  Z^  i is  an 
nxn  zero  matrix. 

The  expansion  of  Eq  (60)  is  straightforward  once 
the  component  matrices  have  been  defined.  The  exponential 
expansion  is  substituted  into  the  sensitivity  transfer 
function  equation.  The  remaining  operations  are  simple 
matrix  algebra. 
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APPENDIX  B 


Minimal  Polynomial 


The  characteristic  polynomial  of  a matrix  M is 


defined  as 


A (X ) = det  (XI-M)  = n (X-X.  ) 

k=l  K 


where  q is  the  number  of  distinct  eigenvalues  of  M and 


n^  is  the  multiplicity  of  X^  in  M.  The  adjoint  matrix  of 


M is  defined  to  be  the  transpose  of  the  matrix  of  cofactors 


of  M (see,  for  example.  Ref  12:211).  The  largest  factor 


that  is  common  to  all  elements  of  the  adjoint  matrix  of 


(XI-M)  is  called  g(X)  and  is  at  least  unity;  that  is. 


b • • • b 
1 • 1 x »n 


adj  (XI-M)  = * 


= g (X)  * 


a a 

n 1 1 • • • n ,n 


The  minimal  polynomial  of  M is  calculated  from 


q m. 

m (x)  - A(x)/g(x)  - n (x-xv) 

k=l  K 


where  m^,  n^  ^m^^l,  is  the  multiplicity  of  X^  in  the 
minimal  polynomial  of  M (Ref  13:65). 
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